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Abstract 


The  Autonomous  Ocean  Sampling  Network  Phase  Two  (AOSN-II)  experiment 
was  conducted  in  and  offshore  from  the  Monterey  Bay  on  the  central  California  coast 
during  July  23  -  September  6,  2003.  The  objective  of  the  experiment  was  to  learn  how  to 
apply  new  tools,  technologies,  and  analysis  techniques  to  adaptively  sample  the  coastal 
ocean  in  a  manner  demonstrably  superior  to  traditional  methodologies,  and  to  use  the 
information  gathered  to  improve  predictive  skill  for  quantities  of  interest  to  end-users. 

The  scientific  goal  was  to  study  the  upwelling/relaxation  cycle  near  an  open  coastal  bay 
in  an  eastern  boundary  current  region,  particularly  as  it  developed  and  spread  from  a 
coastal  headland.  The  suite  of  observational  tools  used  included  a  low-flying  aircraft,  a 
fleet  of  underwater  gliders,  including  several  under  adaptive  autonomous  control,  and 
propeller-driven  AUVs  in  addition  to  moorings,  ships,  and  other  more  traditional 
hardware.  The  data  were  delivered  in  real-time  and  assimilated  into  the  Harvard  Ocean 
Prediction  System  (HOPS),  the  Navy  Coastal  Ocean  Model  (NCOM),  and  the  Jet 
Propulsion  Laboratory  implementation  of  the  Regional  Ocean  Modeling  System  (JPL/ 
ROMS). 

Two  upwelling  events  and  one  relaxation  event  were  sampled  during  the 
experiment.  The  upwelling  in  both  cases  began  when  a  pool  of  cold  water  less  than  13  °C 
appeared  near  Cape  Ano  Nuevo  and  subsequently  spread  offshore  and  southward  across 
the  bay  as  the  equatorward  wind  stress  continued.  The  primary  difference  between  the 
events  was  that  the  first  event  spread  offshore  and  southward,  while  the  second  event 
spread  only  southward  and  not  offshore.  The  difference  is  attributed  to  the  position  and 
strength  of  meanders  and  eddies  of  the  California  Current  System  offshore  which  blocked 
or  steered  the  cold  upwelled  water.  The  space  and  time  scales  of  the  mesoscale  variability 
were  much  shorter  than  have  been  previously  observed  in  deep-water  eddies  offshore. 
Additional  process  studies  are  needed  to  elucidate  the  dynamics  of  the  flow. 
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Introduction 


The  drive  towards  integrated  eoastal-oeean  observing  systems  has  fostered  the 
need  for  basic  research  not  only  on  the  individual  elements  composing  such  systems  but 
also  on  how  to  make  these  elements  function  together  smoothly  as  a  whole.  As  a  step 
towards  providing  improved  coastal  forecasting  over  the  coming  decades,  the  second 
Autonomous  Ocean  Sampling  Network  (AOSN-II)  experiment  was  conducted  in  the 
Monterey  Bay  during  late  July  through  early  September  2003.  The  objective  of  the 
Predictive  Skill  Experiment  was  to  learn  how  to  apply  new  tools,  technologies,  and 
analysis  techniques  to  adaptively  sample  the  coastal  ocean  in  a  manner  demonstrably 
superior  to  traditional  methodologies,  and  to  use  the  information  gathered  to  improve 
predictive  skill  for  quantities  of  interest  to  end-users.  These  quantities  might  include  the 
location  of  fronts  and  boundaries  relevant  to  outfall  dispersal,  particle  trajectories  for 
people  or  pollutants  in  the  water,  and  the  biological  response  to  ocean  physics  in  the  form 
of  bioluminescence  or  harmful  algal  blooms.  The  goal  necessitates  an  underwater, 
through-the-air,  and  over-land  networking  infrastructure  from  which  the  program  takes  its 
name. 


Although  AOSN-II  was  a  novel  collaborative  effort  involving  state-of-the-art 
autonomous  observing  platforms  and  two  data-driven  models  miming  in  real-time  (Chao 
et  al,  this  volume;  Haley  et  al,  2007)  with  data  assimilation  and  ensemble  uncertainty 
predictions  (Lermusiaux,  2007),  the  AOSN-II  project  was  influenced  by  previous  real¬ 
time  applications  of  observing  networks  and  modeling  systems  in  coastal  regions.  Some 
of  the  early  ocean  observing  and  prediction  system  (OOPS)  efforts  are  reviewed  by 
Robinson  et  al  (1998)  and  Robinson  and  Lermusiaux  (2002).  Significant  observing  and 
prediction  coastal  efforts  are  also  described  in  manuscripts  (Lynch  and  Davies,  1995; 
Mooers,  1999;  Pinardi  and  Woods,  2002)  and  review  articles  (e.g.  Robinson  and  Glenn, 
1999;  Dickey,  2003;  Lermusiaux  et  al,  2006).  Without  being  exhaustive,  research  on 
coastal  OOPS  include  scientific  research  exercises  in:  Boston  Harbor  and  Massachusetts 
Bay  (Littoral  Ocean  Observing  and  Prediction  Systems:  Robinson  and  the  LOOPS  group, 
1999;  Lermusiaux,  2001);  in  Georges  Bank  (Lynch  et  al,  2001);  and,  in  the  Martha’s 
Vineyard  region  within  the  Coupled  Marine  Boundary  Layers  and  Air-Sea  Interaction 
Initiative  (CBLAST:  Edson  et  al,  1999;  2006).  AUVs  and  models  have  also  been  jointly 
utilized  within  sustained  efforts  in  the  Gulf  of  Maine  such  as  the  Gulf  of  Maine  Ocean 
Observing  System  (GoMOOS,  Bogden  et  al,  2001)  and  the  Regional  Association  for 
Research  on  the  Gulf  of  Maine  (RARGOM:  Runge,  2005).  Several  of  these  research 
activities  in  the  Northeast  involved  coupled  physical-biological  observing  and  modeling 
research,  from  determining  dominant  scales  to  targeted  process  studies,  e.g.  Ecology  of 
Harmful  Algal  Blooms/Global  Ocean  Ecosystem  Dynamics  (ECOHAB/GLOBEC: 

Anon.,  2002,  2005;  Wiebe  et  al,  2001)  and  LOOPS  (Beskitepe  et  al,  2003),  but  for 
biological-physical  data  collection  and  modeling  in  real-time,  interesting  challenges 
remain  {e.g.  deYoung,  2004).  The  first  Long-Term  Ecosystem  Observatory  (LEO) 
research  with  advanced  ocean  and  observing  prediction  systems  has  been  maintained  off 
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the  coast  of  New  Jersey,  with  numerous  data-model  integrated  coastal  studies  and  results 
(e.g.  Scholfield  et  al,  2002;  Wilkin  et  al,  2005).  Important  research  has  also  being  carried 
out  in  the  Gulf  or  Mexico  and  its  larger  region  of  influence  (e.g.  Mooers  et  al,  2005). 
Numerous  scientific  results  have  also  been  obtained  off  the  coast  of  Oregon  as  a  result  of 
coupled  observation  and  modeling  experiments  (e.g.  Barth  and  Wheeler,  2005).  All  of 
these  research  efforts  are  now  better  coordinated  and  developed  at  the  national  level, 
within  programs  such  as  the  Integrated  and  Sustained  Ocean  Observing  Systems 
(Ocean.US,  2002;  2004).  At  the  international  level,  several  exercises  were  sponsored  in 
multiple  regions  by  the  North  Atlantic  Treaty  Organization  within  the  Rapid 
Environmental  Assessments  and  Rapid  Response  exercises  (e.g.  Bovio  et  al,  2003; 
Robinson  et  al,  2003;  Onken  et  al,  2002,  2005).  Specifically,  the  first  adaptive  sampling 
based  on  real-time  data-driven  uncertainty  predictions  was  carried  out  within  in  a  NATO 
research  exercise  in  the  Strait  of  Sicily  region  (Robinson  et  al,  1999,  Lermusiaux,  1999). 

While  many  of  the  objectives  of  the  AOSN-II  experiment  were  practical  in  nature,  the 
AOSN-II  science  team  formulated  a  list  of  scientific  goals  as  well  including: 

•  Observe  and  predict  the  development  and  movement  of  upwelling  fronts  and 
ocean  eddies 

•  Determine  the  importance  of  the  wind  stress  curl  in  driving  upwelling  and  the 
three-dimensional  circulation  in  the  Monterey  Bay 

•  Study  the  relationship  between  micro-scale  atmospheric  jets  and  the  formation  of 
oceanic  cold  plumes  off  coastal  headlands  and  promontories 

•  Examine  the  nutrient  supply  processes  in  the  Monterey  Bay  including  upwelling, 
advection,  and  mixing 

•  Evaluate  and  improve  numerical  models  of  both  the  atmosphere  and  the  coastal 
ocean. 

The  Monterey  Bay  provides  an  ideal  location  for  such  an  experiment  due  to  its  close 
proximity  to  several  large  research  institutions,  ship  availability,  and  the  wealth  of 
environments  and  scientific  problems  presented.  The  wide  continental  shelves  within  and 
to  the  north  of  the  bay  provide  a  strong  contrast  to  the  nearby  deep  water  within  the 
Monterey  Bay  Submarine  Canyon  (MSC)  and  the  narrower,  steeper  shelf  to  the  south. 

The  two  upwelling  centers  at  Point  Sur  and  Cape  Ano  Nuevo  are  both  within  a  day’s 
steam  by  research  vessel.  Much  of  the  communications  infrastructure  including  undersea 
cables,  semi-permanent  moored  surface  buoys,  and  a  wireless  Internet  network  was 
already  established  by  the  Monterey  Bay  Aquarium  Research  Institute  and  various  local 
Navy  entities.  A  network  of  coastal  HE  radars  was  already  in  place  monitoring  the  surface 
currents  in  real  time.  The  installation  of  the  additional  elements  required  for  the  AOSN-II 
experiment  to  the  existing  observing  system  in  the  bay  was  straightforward. 

The  mesoscale  variability  near  the  Monterey  Bay  can  be  succinctly  described  as 
the  interplay  between  the  upwelling  centers  rooted  at  headlands  to  the  north  and  south  of 
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the  bay  and  a  persistent,  anticyclonic  California  Current  meander  just  offshore  of  the  bay 
itself  (Ramp  et  ah,  2005).  The  upwelling  eenters  are  clearly  driven  by  the  prevailing 
northwesterly  winds  which  are  strongest  during  March  through  June  (Nelson,  1977).  The 
flow  off  Ano  Nuevo  is  sometimes  bifurcated  (Rosenfeld  et  ah,  1994)  with  some  of  the 
cold,  salty  upwelled  water  flowing  offshore  and  some  flowing  south  over  the  Monterey 
Canyon.  This  pattern  seems  to  depend  on  the  wind  stress  duration:  The  cold  water 
spreads  southward  only  after  a  strong,  steady  northwesterly  wind  has  been  blowing  for 
around  6-7  days  (Ramp  et  ah,  2005). 

To  the  south  of  the  bay,  a  strong  upwelling  center  located  at  Point  Sur  generates 
cold  filaments  moving  offshore  in  the  upper  100  m  (Breaker  and  Mooers,  1986;  Tisch  et 
ah,  1992;  Rosenfeld  et  ah,  1994;  Ramp  et  ah,  1997).  On  the  continental  shelf  there, 
currents  were  equatorward  during  northwesterly  winds  and  poleward  during  wind 
reversals  and  relaxations  (Ramp  and  Abbott,  1998).  They  found  that  the  vertical  structure 
of  equatorward  currents  was  consistent  with  the  superposition  of  the  local  wind-forced 
Ekman  spiral  and  the  alongshore  geostrophic  flow  due  to  the  set-down  at  the  coast.  The 
poleward  flows  during  relaxations  were  dynamically  consistent  with  a  poleward 
alongshore  pressure  gradient  force,  as  observed  at  other  locations  off  California  (Winant 
et  ah,  1987;  Lentz,  1987).  These  currents  are  consistent  with  the  view  that  cold  surface 
water  in  the  center  of  the  Monterey  Bay  is  advected  from  the  north  of  the  bay  and  not  the 
south. 


The  anticyclonic  California  Current  (CC)  meander,  also  sometimes  referred  to  as 
the  Monterey  Bay  Eddy  (MBE),  is  a  frequently-observed  feature  of  the  region  during  the 
upwelling  season  (Rosenfeld  et  ah,  1994;  Paduan  and  Rosenfeld,  1996;  Ramp  et  ah, 
2005).  Drifters  suggest  that  this  feature  is  associated  with  the  larger-scale,  meandering 
California  Current  System  and  is  not  locally  generated  (Brink  et  ah,  1991).  Moored 
observations  show  that  the  feature  is  deep,  with  coherent  flows  greater  than  20  cm  s'^ 
exceeding  1000  m  depth  (Ramp  et  ah,  1997).  The  feature  is  warm  and  fresh  relative  to  the 
local  waters,  displaying  the  influence  of  the  Pacific  Subarctic  Water  (Rosenfeld  et  ah, 
1994).  The  MBE  is  thus  an  oceanic  feature  and  is  clearly  not  generated  by  the  local  wind 
stress,  although  it  does  respond  to  it.  Time  series  of  sea  surface  temperature  from  aircraft 
and  AVHRR  images  show  that  the  MBE  moves  rapidly  onshore  during  wind  relaxation 
events  and  retreats  back  offshore  when  the  winds  reintensify  (Rosenfeld  et  ah,  1994; 
Ramp  et  ah,  2005).  Moored  observations  however  have  also  shown  an  onshore  and 
southward  translation  of  the  MBE  without  a  wind  relaxation  (Ramp  et  ah,  1997).  The 
mechanism  for  this  across-shore  eddy  translation  is  thus  not  fully  understood. 

Coursing  poleward  beneath  the  upwelling  centers  and  CC  meander  lies  the 
subsurface  California  Undercurrent  (CUC).  The  CUC  is  observed  over  the  continental 
slope  all  along  the  west  coast  of  the  United  States  (Pierce  et  ah,  2000)  and  is  a  ubiquitous 
feature  of  all  moored  data  sets  both  north  and  south  of  the  Monterey  Bay.  To  the  north 
off  the  Farallon  Islands,  the  undercurrent  was  present  most  of  the  year  but  was  not 
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coherent  with  local  wind  forcing  (Noble  and  Ramp,  2000).  To  the  south  off  Point  Sur,  the 
eurrents  below  100  m  depth  are  more  often  poleward  than  equatorward  and  frequently 
exceed  30  cm  s'^  at  100  m  (Chelton,  1984;  Wickham  et  ah,  1987;  Tisch  et  ah,  1992; 

Ramp  et  ah,  1997).  The  poleward  flow  off  Point  Sur  is  pulse-like  at  very  low  frequencies 
(3-4  months)  (Ramp  et  ah,  1997).  How  this  pulse-like  flow  relates  to  cross-shore 
translations  of  the  MBE  and  potential  “blocking”  of  the  CUC  by  the  eddy  is  of  interest 
but  poorly  understood. 

The  circulation  in  the  bay  itself  is  usually  cyclonic  (Paduan  and  Rosenfeld,  1996; 
Paduan  and  Cook,  1997;  Ramp  et  ah,  2005)  during  both  upwelling  and  downwelling 
eonditions.  During  upwelling  events,  very  warm  (>16“C)  water  is  observed  in  the 
northeast  comer  of  the  bay,  in  the  wind  shadow  behind  the  Santa  Cmz  mountains 
(Graham,  1993;  Graham  and  Largier,  1997;  Ramp  et  ah,  2005).  Model  salinity  fields  at 
the  surface  and  100  m  levels  show  that  the  bay  fills  with  higher-salinity  water  stemming 
from  the  Point  Ano  Nuevo  upwelling  center  under  these  conditions.  During 
downwelling,  the  cyclonic  circulation  was  more  barotropic  over  the  continental  shelf  than 
during  upwelling  events.  The  source  water  for  both  the  surface  and  100  m  levels  was  the 
colder,  fresher  California  Current  water  offshore,  which  had  advected  southward  well 
past  Point  Pinos  during  the  previous  upwelling  event  and  entered  the  bay  from  the  south 
(Ramp  et  ah,  2005) 

Data  and  Methods 

An  impressive  array  of  sampling  platforms  was  assembled  for  AOSN-II  including 
ships,  aircraft,  gliders,  and  propeller-driven  AUVs,  with  each  deployed  to  optimize  its 
sampling  strengths  (Figure  1).  Gliders  are  buoyancy-driven  vehicles  that  are  slow  but  can 
remain  deployed  for  several  months  at  a  time  (Bellingham  et  ah,  this  issue).  Propeller- 
driven  vehicles  can  sustain  speeds  of  6  knots  but  only  for  14  hours  or  less  per 
deployment.  Ships  collect  the  highest-accuracy  data  and  can  collect  water  from 
throughout  the  water  column  for  chemical  analyses.  Aircraft  can  obtain  a  “snapshot”  of 
the  sea  surface  conditions  in  five  hours  or  less,  but  with  less  accuracy  than  ships. 

Five  Spray  gliders  (Sherman  et  ah,  2001;  Davis  et  ah,  2002)  were  deployed 
during  AOSN-II  from  July  23  to  August  24,  2003  along  five  linear  tracks  covering  from 
Point  Ano  Nuevo  in  the  north  to  Point  Sur  in  the  south,  out  to  roughly  100  km  offshore 
(Figure  1).  Beyond  instrument  development  and  testing  operational  capabilities,  the  goal 
for  these  gliders  was  to  contribute  data  for  assimilation  and  synoptic  mapping  and  to 
define  the  scales  and  patterns  of  variability.  These  gliders  collected  a  total  of  2,075 
vertical  profiles  to  400  m  depth  (or  the  bottom)  with  a  nominal  spatial  resolution  of  2.3 
km.  The  resolution  is  a  function  of  bottom  depth  and  higher  resolution  was  obtained  in 
shallower  water.  In  addition  to  temperature,  salinity,  and  water  clarity,  each  glider  also 
returned  an  estimate  of  the  vertically  averaged  horizontal  current  during  each  profile, 
obtained  as  the  difference  between  aetual  GPS  position  and  dead  reckoning.  A  Spray 
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glider  can  typically  make  only  about  22  -  25  cm/s  of  forward  progress  through  the  water, 
which  necessitated  the  across-shore  racetracks  to  avoid  fighting  the  primarily  alongshore 
current  in  the  region.  It  took  each  glider  nominally  4  days  to  transit  one  way  along  each 
100  km  racetrack. 

Overall  Spray  glider  performance  was  excellent  during  AOSN-II.  The  gliders, 
once  deployed,  required  no  servicing  and  lasted  for  the  duration  of  the  experiment  on  one 
set  of  batteries.  The  limiting  factor  on  good  data  collection  turned  out  to  be  biofouling 
rather  than  battery  life.  For  additional  details,  see  the  accompanying  technical  article 
(Bellingham  et  ah,  this  issue). 

Ten  Slocum  gliders  profiled  to  200  m  depth,  and  were  thus  deployed  closer  to 
shore  than  the  Sprays  (Figure  1).  Four  Slocum  gliders  ran  closed  boxes  (or  “cells”)  rather 
than  linear  tracks,  and  were  equipped  with  fluorescence  and  PAR  (photosynthetically 
available  radiation)  sensors  in  addition  to  the  variables  sensed  by  the  Spray.  A  total  of 
11,341  profiles  yielding  nominally  800  m  along-track  resolution  were  obtained  over  the 
life  of  the  experiment.  The  endurance  of  the  Slocums  was  11-14  days,  about  half  that 
expected,  and  some  recovery  and  re-deployment  was  thus  required  to  cover  the  lifetime 
of  the  experiment.  The  Slocums  also  experience  problems  with  biofouling,  but  the  shorter 
deployments  allowed  for  cleaning  of  the  sensors  to  minimize  this  problem. 

One  of  the  most  innovative  aspects  of  the  AOSN-II  experiment  was  adaptive 
sampling  of  ocean  features  using  gliders  (Leonard  et  ah,  2007).  Three  to  four  Slocum 
gliders  were  used  for  this  purpose.  Three  demonstrations  of  autonomous  adaptive  control 
were  successfully  completed,  one  in  the  northern  part  of  the  bay  across  an  upwelling 
front,  the  second  in  the  southern  part  of  the  bay,  and  the  third  following  a  surface  drifter 
(Fiorelli  et  ah,  2006).  For  the  first  experiment  during  August  6-7,  2003,  three  gliders 
maintained  a  triangular  pattern  3  km  on  a  side  while  propagating  NW  along  a  linear  path. 
The  triangular  formation  was  chosen  so  that  in-situ  estimates  of  gradients  (e.g.,  in 
temperature)  along  the  path  of  the  center  of  the  formation  could  be  computed  in  near-real 
time  from  the  gliders’  scalar  (temperature)  measurements.  The  gliders  could  be 
programmed  to  follow  the  gradient  if  desired.  The  3  km  sensing  array  resolution  was 
designed  with  respect  to  previously  computed  scales  in  Monterey  Bay.  The  formation 
error  was  423  m,  or  14%  of  the  desired  spacing  with  a  standard  deviation  of  159  m.  For 
the  second  demonstration  on  August  16-17,  2003,  some  additional  difficulty  was  added. 
The  triangle  formation  was  asked  to  travel  a  zigzag  path  across  an  upwelling  front,  and  to 
shrink  the  size  of  the  formation  from  6  to  3  km  on  a  side  enroute.  The  objective  was  to 
demonstrate  the  ability  of  the  glider  control  system  to  adjust  the  formation  size,  i.e.,  the 
sensing  array  resolution,  in  response  to  changing  scales  in  the  sampled  field(s)  and  to 
provide  gradient  estimates  at  different  resolutions.  This  exercise  was  also  a  success:  The 
glider  formation  error  was  a  function  of  the  size  of  the  triangle,  about  394  m  for  the  6  km 
triangle  and  270  m  for  the  3  km  one.  Finally,  on  August  23,  2003,  a  single  glider  was 
asked  to  follow  a  surface  drifter  in  real  time,  making  zigzags  across  its  projected  path. 
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The  goal  was  to  demonstrate  that  a  glider  or  glider  formation  eould  colleet  sealar 
measurements  (or  gradient  estimates)  both  aeross  and  below  traeer  paths.  This  demo  was 
qualitatively  a  sueeess,  but  there  were  some  problems  based  on  the  different  veloeity 
fields  sensed  by  the  two  instruments:  The  drifter  measures  the  surfaee  veloeity  while  the 
glider  measures  the  current  vertically  integrated  over  the  up/down  path,  making 
coincident  path  prediction  difficult.  There  are  other  (easier)  ways  to  do  this,  such  as 
having  the  glider  follow  the  drifter  with  a  two-hour  lag,  which  may  be  more  useful  in 
different  applications. 

Propeller-driven  AUVs  such  as  the  commercially  available  REMUS,  the  MBARI 
manufactured  DORADO,  and  the  NPS  in-house  vehicle  ARIES  typically  have  about  14 
hours  duration  at  4-6  knots.  The  vehicles  also  have  a  much  larger  payload  than  gliders 
and  can  thus  carry  more  sophisticated  sensors.  The  largest  number  of  transects  were  run 
by  the  California  Polytechnic  State  University  (CalPoly)  REMUS  vehicle  along  two 
transects  in  the  northern  half  of  the  Monterey  Bay.  The  goal  of  sensing  bioluminescence 
dictated  that  these  REMUS  runs  be  done  in  the  dark,  and  therefore  nine  seven-hour 
missions  were  completed  along  45  km  transects  between  9  pm  and  4  am  each  night  (blue 
“V”  in  Figure  1).  These  operations  began  near  Santa  Cruz,  ran  out  to  MBARI  surface 
buoy  Ml,  then  returned  to  the  NW  comer  of  the  bay  (Moline  et  ah,  2005;  Shulman  et  ah, 
2005).  The  vehicle  followed  a  saw-tooth  pattern  to  40  m  depth  sampling  with  a  CTD, 
transmissometer,  fluorometer,  and  ADCP  in  addition  to  the  bioluminescence.  Four 
transects  were  also  run  from  Monterey  to  Ml  using  the  NPS  REMUS  vehicle,  sampling 
temperature  only.  The  NPS  ARIES  vehicle  was  used  exclusively  to  fetch  data  from  a 
bottom-mounted  ADCP  using  acoustic  modems  and  did  not  collect  environmental  data. 

The  largest  AUV  used  during  AOSN-II  was  the  MBARI  DORADO  AUV.  This 
vehicle  ran  along  a  line  between  Moss  Eanding  and  MBARI  buoy  M2,  profiling  to  250  m 
depth  (Figure  1).  Instmments  on  board  included  a  CTD,  fluorometer,  oxygen  and  nitrate 
sensors,  bioluminescence,  and  ADCP.  The  mn  was  made  operationally  every  three  weeks 
but  the  frequency  was  increased  during  August  2003  to  1  section  per  week.  Two 
additional  DORADO  transects  were  also  conducted  offshore.  A  clear  advantage  of  this 
vehicle  was  the  sampling  depth  which  was  very  helpful  for  model  data  assimilation.  Its 
large  size  requires  a  specialized  vehicle  for  recovery. 

The  airborne  measurements  during  AOSN-II  were  obtained  using  the  TWIN 
OTTER  aircraft  owned  and  operated  by  the  Center  for  Interdisciplinary  Remotely-Piloted 
Aircraft  Studies  (CIRPAS)  at  the  Naval  Postgraduate  School.  Fourteen  flights  were 
conducted  between  August  4  and  September  5,  2003  at  an  altitude  of  nominally  33  m 
above  the  sea  surface.  The  low  altitude  was  chosen  to  pass  beneath  the  stratus  cloud  deck 
that  often  covers  the  Monterey  Bay  in  summer,  and  to  minimize  the  vertical  divergence  in 
the  heat  and  momentum  flux  observations.  The  flight  path  (Figure  1,  blue  line)  covered 
all  of  the  bay  and  some  distance  to  the  north  and  south  out  to  100  km  offshore.  Two 
vertical  saw-tooth  paths  to  1500  m  were  also  conducted  along  two  slightly  longer 
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transects  at  the  northern  end  of  the  flight  path  and  along  the  center  of  the  bay  (Figure  1). 
These  paths  were  to  sample  the  three-dimensional  structure  of  the  marine  atmospheric 
boundary  layer  over  the  bay  (Kalogiros  et  ah,  2006;  Wang  et  ah,  2007).  The  atmospheric 
parameters  sampled  by  the  aircraft  included  air  temperature  and  dew  point,  atmospheric 
pressure,  wind  speed  and  direction,  the  turbulent  fluxes  of  heat  and  momentum,  and  total 
particle  number,  as  well  as  aerosol  and  cloud/fog  droplet  size  distributions.  In  the  ocean, 
sea  surface  temperature  (SST),  ocean-leaving  radiance  at  193  wavelengths,  sun  glint,  and 
surface  roughness  were  observed.  The  SST,  air  temperature  and  dew  point,  and  wind  data 
were  made  available  for  model  assimilation  in  near-real  time,  while  the  other  parameters 
required  additional  post-processing. 

Coastal  HF  radars  for  mapping  ocean  surface  currents  have  been  proliferating 
along  the  central  coast  for  a  period  of  years,  starting  with  three  standard-range  (13  MHz) 
systems  in  the  Monterey  Bay  located  at  Point  Pinos,  Moss  Landing,  and  Santa  Cruz.  The 
network  has  recently  been  extended  to  the  north  with  long-range  (5  MHz)  systems  at  Big 
Creek,  Pescadero,  and  Montara.  All  the  bay  systems  were  operational  during  AOSN-II 
and  provided  surface  currents  in  2-3  km  bins  for  model  data  assimilation  (Paduan  and 
Shulman,  2004)  and  Lagrangian  Coherent  Structures  (LCS)  computations. 

The  Numerical  Models 

A  key  aspect  of  AOSN-II  program  was  to  assimilate  the  data  collected  in  near- 
real  time  into  mesoscale  ocean  models  to  forecast  conditions  for  the  following  day,  and 
adapt  the  sampling  scheme  accordingly.  In  addition  to  real  time  operations,  the  models 
were  also  used  to  help  elucidate  the  dynamics  of  the  flow.  Three  models  were  used  for 
these  purposes,  namely  the  Harvard  Ocean  Prediction  System  (HOPS)  (Robinson,  1999), 
the  Jet  Propulsion  Laboratory  implementation  of  the  Regional  Ocean  Modeling  System 
(JPL/ROMS)  (Shchepetkin  and  McWilliams,  2004),  and  the  Navy  Coastal  Ocean  Model  / 
Innovative  Coastal  Ocean  Observing  Network  (NCOM/ICON)  (Shulman  et  ah,  2002; 
2004).  While  conceptually  similar  in  scale  (1.5  -  2.5  km  in  the  horizontal  x  32  sigma 
levels  in  the  vertical)  and  operation,  the  models  differed  significantly  in  terms  of  the 
mixing  schemes  employed,  initial  conditions  used,  spin-up  time,  and  open  boundary 
conditions  provided  by  the  larger-scale  models  they  were  nested  within.  All  three  models 
were  forced  using  the  Coupled  Air  Ocean  Modeling  and  Prediction  System  (COAMPS™) 
surface  wind  stresses  and  heat  fluxes,  and  assimilated  the  aircraft  SST,  glider  T  and  S, 
and  AUV  data.  The  JPL/ROMS  model  additionally  assimilated  T  and  S  data  from 
MBARI  moorings  Ml  and  M  2.  None  of  the  models  assimilated  velocity  from  moorings 
or  HF  radar.  A  table  with  additional  details  of  the  three  models,  the  data  assimilation 
schemes  used,  and  how  they  were  forced  may  be  found  on  the  web  at:  (http:// 
www.princeton.edu/~dcsl/asap/ 

ASAP_ROMS_HOPS_NCOMICON_summary_MB06.doc). 

Results 
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Atmospheric  Forcing 


The  weather  presented  an  interesting  mix  of  eonditions  whieh  were  well  suited  to 
the  experimental  objeetives  during  August  2003.  Time  series  from  MBARI  buoy  M2, 
loeated  at  36°  42'  N,  122°  23'  W  at  the  mouth  of  the  bay,  show  that  the  first  two  flights 
during  August  4-5  took  plaee  during  a  wind  relaxation  event  when  winds  were  weakly 
poleward  at  less  than  5  m  s'^  (Figure  2).  An  extended  period  of  strong  (eommonly  10  m 
s'^),  steady,  upwelling-favorable  northwesterly  winds  followed  during  August  6-19, 
which  were  quite  typical  for  this  time  of  year.  Five  TWIN  OTTER  flights  sampled  this 
first  upwelling  event.  Another  poleward  event  took  place  during  August  20-21,  followed 
by  a  second  upwelling  event  from  August  23-31.  Three  flights  were  executed  during  the 
relaxation  and  two  during  the  second  upwelling  event.  September  1-3  was  very  calm 
and  warm  in  the  Monterey  Bay  and  was  characterized  by  nearly  fiat  seas  and  strong 
surface  heating.  This  time  period  was  unfortunately  not  sampled  by  the  aircraft.  During 
September  4-6,  weakly  upwelling-favorable  winds  returned  to  the  area  once  again,  and 
the  final  two  aircraft  over-flights  were  conducted. 

Since  the  spatial  scales  of  the  atmospheric  forcing  are  large  compared  to  the  size 
of  the  Monterey  Bay  (22  x  40  km),  the  winds  observed  at  buoy  M2  are  usually 
representative  of  the  offshore  winds  elsewhere  along  the  coast  (Figure  3).  Within  the  bay 
itself,  however,  the  surrounding  mountains  create  high  spatial  variability  in  the  surface 
wind  stress.  Of  particular  interest  are  the  weak  winds  behind  the  mountains  in  the  NE 
corner  of  the  bay,  and  a  small-scale  atmospheric  jet  coming  through  a  gap  in  the 
mountains  off  Santa  Cruz  (Figure  3).  This  figure  from  August  15  was  sampled  nine  days 
into  the  upwelling  event  and  was  typical  of  the  others  during  this  time.  The  aircraft  and 
the  Coupled  Ocean- Atmosphere  Mesoscale  Prediction  System  (CO AMPS)  (Doyle  et  ah, 
this  issue)  both  show  the  upwelling  favorable  winds  turning  slightly  into  the  bay.  Both 
also  show  the  weak  (blue)  winds  in  the  NE  comer  of  the  bay.  The  wind  jet  off  Santa  Cmz 
was  slightly  clearer  in  the  model  output  than  the  observational  data,  as  was  the  strong 
positive  curl  on  the  inshore  side  of  the  jet  (Figure  3,  bottom).  The  location  and  scale  of 
this  feature  closely  resembled  a  similar  one  observed  on  August  17,  2000  (Ramp  et  ah, 
2005).  While  that  jet  was  noticeably  hot  and  dry  (their  Figure  4)  this  was  less  obvious  in 
the  2003  data.  This  is  perhaps  due  to  the  higher  altitude  (133  m)  of  the  year  2000  flights. 
The  2003  atmospheric  temperatures  at  33  m  (not  shown)  more  closely  followed  the  sea 
surface  temperature  and  were  cool  rather  than  warm.  There  was  a  minimum  in  the  dew 
point  temperature  just  off  Santa  Cmz  however,  indicating  slightly  dryer  conditions  there. 
A  typical  poleward  wind  event  from  August  20,  2003  (Figure  4)  showed  overall  weaker 
winds,  weaker  gradients,  and  less  wind  stress  curl  than  a  typical  equatorward  wind  event. 
The  wind  stress  was  still  weaker  in  the  NE  corner  of  the  bay  than  elsewhere,  and  the  wind 
turned  into  the  bay  from  the  south  rather  than  the  north.  The  model/data  comparisons  are 
discussed  in  greater  detail  elsewhere  (Paduan  et  ah,  this  issue). 


11 


Development  and  relaxation  of  the  upwelling  plume 


A  time-series  of  twelve  realizations  of  SST  and  wind  as  sampled  by  the  aircraft 
shows  the  ocean’s  mesoscale  response  to  the  equatorward  and  poleward  wind  stress 
events  during  summer  2003  (Figure  5a-l).  August  4  was  the  first  day  of  a  relaxation  event 
following  a  weak,  2-day  upwelling  event.  The  winds  were  weak  and  disorganized  and 
some  traces  of  coastal  upwelling  were  still  visible  (<14  °C)  off  Point  Ano  Nuevo.  The 
warmest  water  exceeded  17  °C  offshore  and  in  the  southern  portion  of  the  Monterey  Bay. 
By  the  following  afternoon  (Figure  5b),  the  entire  sea  surface  had  warmed  by  about  2  °C 
under  weak,  slightly  onshore  winds,  with  the  exception  of  off  Point  Sur  where  the  SST 
remained  about  the  same.  There  was  some  indication  of  locally  upwelling  favorable 
winds  off  Point  Sur  which  perhaps  balanced  the  local  surface  heating. 

The  next  five  panels  (Figure  5c-g)  show  the  development  of  the  largest  upwelling 
event.  The  August  6  image  was  obtained  on  the  first  day  of  strong  (10  m  s'^)  upwelling- 
favorable  winds.  The  upwelling  was  just  starting  (yellow  contour  <15  °C)  at  both  Point 
Ano  Nuevo  and  Point  Sur.  Note  that  the  aircraft  survey  was  designed  to  focus  on  the  Ano 
Nuevo  region  and  most  of  the  Point  Sur  upwelling  plume  was  located  to  the  south  of  the 
flight  path  and  was  not  observed.  The  discussion  hereafter  focuses  on  the  growth  and 
spreading  of  the  Ano  Nuevo  plume.  Four  days  later  (Figure  5d)  the  upwelling  center  had 
enlarged  and  cooled  considerably  to  less  than  13  °C.  There  was  some  evidence  of  plume 
bifurcation  (Rosenfeld  et  ah,  1994)  with  some  of  the  14  degree  water  heading  offshore 
and  some  heading  south.  By  the  next  day  however  (Figure  5e),  southward  advection 
seemed  to  be  winning  with  all  the  13  and  14  °C  water  heading  south  and  connecting  up 
with  the  cold  water  from  the  Point  Sur  upwelling  center.  This  lead  to  a  detached  parcel  of 
13  °C  water  which  advected  south  on  August  13  (Figure  5f).  Note  that  the  CC  meander, 
as  indicated  by  the  17  °C  contour,  remained  well  offshore  during  this  entire  event,  which 
allowed  the  cold  water  to  follow  an  unobstructed  path  towards  the  south.  On  August  15, 
following  nine  days  of  continuous  upwelling,  the  winds  were  blowing  at  15  m  s' ^  and  the 
coldest  water  (<12  °C)  was  observed  off  both  Ano  Nuevo  and  Point  Sur.  The  cold  water 
from  the  two  centers  merged  together  and  formed  a  continuous  band  of  minimum  SST 
roughly  along  122°  06'W  that  separated  the  warm  water  offshore  from  the  also  warm 
water  in  the  Monterey  Bay. 

The  two-week  upwelling  event  was  followed  by  three  days  of  poleward  winds  at 
about  8  m  s''on  August  20-22  (Figure  2,  Figures  5h-j).  The  ocean’s  adjustment  to  the 
change  in  surface  forcing  was  remarkably  rapid  and  shows  two  noteworthy  features 
(Figure  5h).  First,  the  boundary  between  the  warm  California  Current  water  moved 
onshore  to  within  just  a  few  kilometers  of  the  coast  near  Davenport,  just  north  of  Santa 
Cruz,  CA.  This  dramatically  reduced  the  amount  of  cold  water  at  the  surface  in  the 
northern  portion  of  the  domain.  Second,  a  cold  “squirt”  formed  in  the  southern  half  of  the 
domain,  centered  near  Point  Pinos  in  Monterey.  Since  the  wind  generally  reverses  from 
south  to  north  along  the  coast,  we  speculate  that  this  feature  resulted  from  convergence 
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between  poleward  currents  along  the  southern  Big  Sur  coast  and  currents  which  were  still 
equatorward  to  the  north  of  this  region.  This  feature  was  still  observable  as  a  lens  of  15 
°C  water  55  km  offshore  on  August  2P‘  (Figure  5i)  but  was  ultimately  pushed  back 
onshore  by  the  continuing  incursion  of  CC  water  on  the  22"^^  (Figure  5j). 

The  final  two  panels  of  Figure  5  show  the  return  to  upwelling  conditions,  as 
represented  by  two  realizations  on  August  25  and  29  (Figure  5k,  1).  The  25*'^  was  just  two 
days  into  the  event  and  not  surprisingly  resembles  August  6*'^,  with  weak  upwelling  and 
<15  °C  water  off  Cape  Ano  Nuevo.  August  29,  six  days  into  the  event  with  winds  >15  m 
s'^  closely  resembles  August  15*''.  One  significant  difference  would  be  that  the  CC  was 
closer  to  shore  on  the  29**'  than  the  15*'',  which  confined  the  cold  water  closer  to  shore  and 
resulted  in  strictly  southward  rather  than  offshore  spreading.  These  differences  in  the  two 
upwelling  events  were  also  evident  in  the  output  from  the  NCOM/ICON  model  (Figure 
6).  The  NCOM-ICON  model  is  based  on  the  Navy  Coastal  Ocean  Model,  and  is  triply 
nested  inside  of  the  global  and  regional  (California  Current)  NCOM-based  models 
(Shulman  et  ah,  2004).  The  model  is  called  NCOM  ICON  due  to  the  fact  that  initial 
development  of  the  model  started  under  the  National  Oceanic  Partnership  Program 
(NOPP)  Innovative  Coastal-Ocean  Observing  Network  (ICON)  project  (Ramp  et  ah, 
2005,  Shulman  et  ah,  2005).  The  top  and  bottom  panels  of  Figure  6  are  a  remarkably 
good  match  for  figures  5g  and  51  and  show  a  more  pronounced  plume  of  cold,  salty  water 
moving  offshore  off  Ano  Nuevo  during  the  August  15"'  event  than  during  the  August  29**' 
event.  The  model  fields  show  the  Monterey  Bay  eddy  to  be  much  closer  to  shore  during 
the  latter  event,  centered  near  36.6°N,  122. 5°W  with  poleward  flow  on  the  offshore  side. 
These  eddy  circulations  were  responsible  for  the  offshore  advection  of  the  cold,  salty 
plumes  off  Point  Pinos  and  Point  Sur,  and  for  suppressing  the  offshore  flow  off  Point  Ano 
Nuevo  on  August  29. 

While  the  previous  discussion  focuses  on  the  evolution  of  sea  surface 
temperature,  some  insight  into  the  vertical  structure  of  the  evolving  upwelling  plume  was 
also  obtained  from  the  glider  formation  flying  experiments  (GFFEs)  described  earlier. 

The  first  GFFE  on  August  6-7  sampled  the  front  when  it  was  young  and  close  to  the 
source  at  Point  Ano  Nuevo  and  the  second  GEFE  on  August  1 6  sampled  the  front  when  it 
was  older  and  had  advected  further  to  the  south  across  the  mouth  of  the  Monterey  Bay.  In 
response  to  the  strong  upwelling-favorable  winds  measured  on  August  6  (Figure  5c), 
three  gliders  were  directed  that  same  day  to  travel  northwest,  in  an  equilateral  triangle 
with  3  km  length  sides,  towards  the  Point  Ano  Nuevo  region  to  sample  along  the 
southwest  side  of  the  upwelling  front  (Figure  7).  Note  that  this  formation  flew  against  the 
prevailing  current  as  sampled  by  all  the  gliders  (blue  arrows  in  the  figure),  which  would 
not  be  possible  using  conventional  profiling  floats  or  drifters.  The  temperature  gradient 
estimate  along  the  path  provides  a  higher  resolution  complement  to  the  SST  (Figure  5c) 
and  in  the  corresponding  AVHRR  plots  (not  shown).  The  temperature  gradients  had  the 
same  sign  at  both  depths  shown  (5  m  and  30  m),  with  warmer  water  offshore  and 
approximately  1  °C  cooler  water  inshore.  The  gradient  was  slightly  stronger  at  30  m. 


13 


below  the  surface  mixed  layer  where  spatial  homogeneity  tends  to  weaken  the  horizontal 
gradients. 


During  the  second  GFFE  on  August  16-17  the  gliders  traveled  from  near  the 
center  of  the  bay  towards  the  southeast  for  about  12  hours  and  then  turned  east  and  then 
northeast  for  a  few  more  hours  (Figure  8).  In  this  case  the  gliders  moved  mostly  with  the 
prevailing  current  (blue  arrows)  but  also  across  and  slightly  against  the  flow  towards  the 
end.  The  temperature  gradients  in  the  horizontal  plane  were  calculated  from  the  glider- 
observed  temperatures  at  5  and  30  m  depth  along  the  path  as  in  the  August  6  experiment. 
During  the  first  12  hours  the  glider  triangle  was  straddling  the  in-shore  boundary  of  the 
upwelling  front,  i.e.,  the  band  of  minimum  SST  near  122°  06'W  (Figure  5g).  At  5  m 
depth,  the  glider(s)  furthest  offshore  measured  cool  temperatures  of  1 1.5  to  12  °C  in  the 
extended  plume  while  the  gliders  to  the  northwest  measured  warmer  temperatures  around 
13  °C  (Figure  8).  At  30  m,  the  temperature  was  nearly  uniform  (10  °C,  dark  blue)  over 
much  of  the  glider  formation  path,  but  with  a  weak  gradient  towards  the  northwest, 
opposite  the  surface  layer.  This  indicates  that  the  upwelling  filament  was  quite  shallow, 
less  than  30  m  deep,  after  it  had  advected  southward  across  the  Monterey  Bay.  By  45  m 
depth  (not  shown)  the  temperature  gradient  was  very  small. 


Offshore  Structure 

A  central  AOSN-II  issue  is  how  well  the  observing  array  can  describe  mesoscale 
ocean  structure  over  the  region  of  interest.  Data  assimilating  models  can  be  powerful 
tools  for  synthesizing  the  scattered  and  diverse  AOSN-II  observations  into  a  coherent 
picture,  but  these  models  can  only  describe  that  range  of  phenomena  that  approximately 
obey  the  model  equations  at  the  resolution  they  are  solved.  Little  is  known  about  what 
observations  are  needed  to  support  an  accurate  and  coherent  synthesis  of  mesoscale 
coastal  structure  or  how  these  data  should  be  statistically  characterized  for  the 
assimilation  process.  The  observing  power  of  the  AOSN-II  array  is  described  by  data 
collected  over  two  separate  one-week  periods  under  contrasting  wind  forcing  (Figures  9  - 
1 1).  In  the  first  week  the  offshore  sampling  density  was  not  yet  very  uniform  since  the 
gliders  were  just  getting  on  their  assigned  tracks  (Figure  9a).  During  the  second  week 
some  gliders  were  moved  off  their  assigned  tracks  by  currents  (Figure  9b). 

Winds  over  the  last  week  of  July  and  through  August  6**^  were  weak  and  variable 
(Figure  2).  While  the  moorings  provide  time  series,  the  other  platforms  provide 
instantaneous  samples  along  tracks  and  it  takes  about  a  week  for  them  to  map  out  the 
region.  Even  without  upwelling-favorable  winds  for  a  week  or  more,  near-surface 
temperatures  are  2°C  or  more  colder  near  the  coast  than  50-100  km  offshore.  The  mean 
alongshore-  and  time-averaged  across-shore  temperature  section  for  this  time  period 
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(Figure  10a)  shows  the  near-surface  temperature  gradient  in  Figure  9a  to  be  associated 
with  isotherm  slopes  (and  isopycnal  slopes  not  shown)  in  the  upper  50-100  m  that  support 
geostrophic  shear  with  equatorward  currents  increasing  near  the  surface.  Below  this,  the 
isotherm  and  isopycnal  slopes  are  reversed,  producing  geostrophic  shear  with  poleward 
flow  increasing  above  100  m. 

Upwelling  favorable  winds  began  August  6  and  persisted  through  August  19 
(Figure  2).  Compared  with  the  preceding  calm  period,  near-surface  temperature  cooled 
about  2°C  within  10  km  of  the  coast  as  the  surface  layer  was  blown  offshore  (Figure  9b). 
This  is  also  evident  in  the  up-warping  of  the  isotherms  within  40  km  of  the  coast  (Figure 
10b).  In  contrast,  the  temperature  section  shows  that  at  distances  of  50  km  or  more  from 
shore  the  near-surface  temperature  increased  up  to  1°C  as  a  warm  surface  layer  seems  to 
have  accumulated  in  a  weakly  stratified  deep  mixed  layer  (Figure  10b).  Below  100  m,  no 
large-scale  change  in  temperature  or  temperature  gradient  can  be  detected.  The 
accumulation  of  warm  water  offshore  under  upwelling-favorable  winds  is  surprising 
since  there  is  no  evidence  that  the  wind  strength  or,  consequently,  the  Ekman  transport 
decreases  offshore  or  that  surface  warming  from  above  increases. 

The  offshore  Spray  gliders  sample  depth-averaged  currents  to  400  m  and  show  a 
general,  if  structured,  poleward  flow  offshore  in  both  time  periods  (Figure  9).  This 
manifestation  of  the  California  Undercurrent  was  strongest  during  the  weak-wind  period 
when  it  took  the  form  of  a  20-km  wide  jet  that  meandered  from  25  km  offshore  crossing 
the  southernmost  glider  line  to  65  km  offshore  on  the  northern  line  where  it  exhibited  the 
strongest  currents  of  nearly  30  cm/s  (Figure  9a).  Equatorward  wind  in  the  second  period 
significantly  reduced  the  poleward  signatures  in  depth-averaged  velocity  (Figure  9b). 
Using  depth-  and  alongshore-averaged  current  as  a  function  of  offshore  distance  to 
reference  the  alongshore  geostrophic  shear  computed  from  alongshore-  and  time- 
averaged  sections  of  temperature  (Figure  10)  and  salinity  allows  cross-sections  of 
absolute  alongshore  velocity  to  be  produced.  These  sections  show  that  the  decrease  in 
poleward  depth-average  current  from  the  first  plotted  week  to  the  second  resulted  from 
offshore  changes  in  geostrophic  shear  of  alongshore  current  in  the  upper  150  m  (Figure 
11).  Also  seen  in  Figure  1 1  is  how  the  equatorward  upwelling  jet  in  the  second  time 
period  is  largely  confined  to  the  20  km  closest  to  shore.  More  surprisingly,  aside  from 
these  features,  the  general  pattern  of  alongshore  flow  is  moderately  complex  but 
remarkably  similar  between  the  two  time  periods  even  though  no  observation  contributes 
to  both  time  averages.  It  remains  to  be  seen  if  this  structure  persisted  for  more  than  16 
days. 


Beyond  the  alongshore-  and  time-averaged  patterns  described  in  Figures  10  and 
11,  Figure  9  shows  complex  and  energetic  smaller  scale  features  in  both  near-surface 
temperature  and  depth-integrated  currents.  Some  of  the  complexity  derives  from  plotting 
a  full  week  of  data  together  while  other  features  may  be  a  consequence  of  internal  tides 
that  are  substantially  aliased  by  shipboard  and  glider  sampling.  But  some  of  the  smaller 
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scale  features  have  the  appearance  of  eddies  and  fronts.  Notable  are  oppositely  rotating 
eddies  found  near  the  mouth  of  Monterey  Bay.  During  the  first  week  (Figure  9a)  the  eddy 
is  eyclonic  with  a  elear  poleward  flow  along  the  coast  of  Monterey  Bay  while  it  is  anti- 
cyclonic  in  the  August  8-15  upwelling  period  (Figure  9b).  The  anti-cyclonic  eddy  is 
implicated  as  a  cause  of  the  spread  of  eold  water  away  from  Ft.  Pinos  (the  southern  end 
of  the  Bay  mouth)  seen  in  aireraft  SST  as  shown  in  Figures  5f-h  spanning  August  13-20. 

The  ability  of  the  AOSN-II  sampling  array  to  describe  the  smaller  scale  features 
hinted  at  in  Figure  9  depends  on  the  spatial  and  temporal  seales  of  those  features. 

Roughly  put,  objeetive  mapping  (Bretherton  et  ah,  1975)  shows  that  mapping  aeouraey  in 
the  absence  of  assimilation  into  a  dynamical  model  depends  on  the  number  of 
observations  within  sampling  volumes  whose  dimension  is  set  by  the  field’s  space  and 
time  seales.  If  there  are  enough  points  in  eaeh  sueh  volume  to  suppress  small-seale  noise 
the  field  ean  be  mapped.  In  plarming  the  AOSN-II  field  trial  the  seales  of  coastal 
variability  were  not  well  known.  Figure  12  shows  the  temporally  lagged  correlation  of 
nearly  co-located  temperature  measurements  and  the  horizontally  isotropic  spatially 
lagged  eorrelation  of  nearly  simultaneous  measurements  at  the  same  depth  based  on  all 
glider,  shipboard  and  mooring  data  from  AOSN-II. 

Several  aspects  of  the  measured  scale  structure  deserve  comment.  First,  the  spatial 
and  temporal  1/e-folding  scales  of  15  km  and  2  days,  respeetively,  are  nearly  a  factor  of 
10  smaller  than  the  seales  of  typical  deep-water  mesoscale  eddies,  which  are  commonly 
100  km  and  30  days.  Second,  the  correlation  values  appear  to  inerease  rapidly  as  both  the 
spatial  and  temporal  separations  become  small.  Indeed,  at  the  smallest  separations  of  0.5 
days  and  5  km  the  correlation  is  about  0.8,  suggesting  that  20%  of  the  observed 
temperature  variance  has  smaller  scales  than  this.  Third,  the  correlation  shapes  are 
approximately  the  same  for  all  depths  shown  (and  all  others  that  are  not  shown).  A  more 
complete  discussion  of  the  statistics  of  the  AOSN-II  variability  is  given  by  Davis  et  al. 
(2007)  but  the  main  point  is  that  the  results  in  Figure  12  show  that  the  sampling  array  for 
the  five  offshore  Spray  gliders  was  too  thin  to  map  any  but  the  largest-seale  eomponents 
of  the  observed  fields. 

Multi-scale  regional  dynamics 

The  multi-seale  dynamieal  interaction  between  the  large-seale  and  mesoscale 
variabilities  of  the  August  2003  cireulation  and  their  relationship  with  the  upwelling- 
related  events  have  been  investigated  using  the  methodology  multi-scale  energy  and 
vorticity  analysis  (MS-EVA),  and  the  MS-EVA-based  localized  instability  theory 
developed  by  Eiang  and  Robinson  (2005a,  2005b).  MS-EVA  is  built  on  the  basis  of  the 
multi-scale  window  transform,  a  recently  established  mathematieal  apparatus  by  Eiang 
and  Anderson  (2007).  In  this  framework,  the  multi-seale  processes  in  a  given  funetional 
space  are  represented  on  its  scale  windows,  or  mutually  orthogonal  subspaces,  each  with 
an  exclusive  range  of  time  scales.  The  objeetive  of  the  MS-EVA  is,  accordingly,  to 
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investigate  how  energy  and  enstrophy  are  transported  within  their  individual  scale 
windows,  redistributed  through  the  interaction  between  the  windows,  and  converted  and 
dissipated  when  a  system  is  steered  forward.  The  redistribution  is  characterized  by  a 
concept  called  perfect  transfer  which  has  been  rigorously  coimected  to  the  two  important 
processes,  barotropic  instability  and  baroclinic  instability,  in  geophysical  fluid  dynamics. 

The  AOSN-II  Harvard  Ocean  Prediction  System  (HOPS)  model  fields  have  been 
analyzed  in  three  mutually  exclusive  time  sub-spaces:  a  large-scale  window  (>  8  days),  a 
mesoscale  window  (0.5  -  8  days),  and  a  sub-mesoscale  window  (<  0.5  days).  The 
analysis  has  been  performed  during  the  upwelling-relaxation-upwelling  cycle  of  winds 
shown  in  Figure  2.  Upwelling  favorable  winds  occur  from  6-19  August  and  23-31  August 
and  the  relaxation  period  is  19-22  August.  The  dynamics  of  the  fields  are  driven 
externally  by  wind  and  surface  fluxes  and  are  characterized  internally  by  the  transfers  of 
available  potential  energy  (APE)  and  kinetic  energy  (KE)  from  the  large  scale  to  the 
mesoscale  both  (i)  during  upwelling  and  (ii)  during  relaxation.  Off  of  Point  Sur  the  wind 
destabilizes  the  system  directly  during  upwelling,  with  wind  energy  captured  within  the 
large-scale  window  and  continually  released  to  the  mesoscale.  In  contrast,  near  and  in 
Monterey  Bay,  the  wind  tends  to  stabilize  the  southward  coastal  current,  energy  is  stored 
in  the  large-scale,  and  instability  occurs  only  when  the  external  wind  constraint  is 
relaxed. 

Figure  13  illustrates  the  available  potential  and  kinetic  energy  transfer  from  the 
large-scale  window  to  the  mesoscale  window  at  four  day  intervals  from  13-25  August  and 
Figure  14  schematizes  the  energy  transfer  process.  During  (i)  upwelling  (13,  17  and  25 
August)  the  greatest  direct  transfer  of  APE  and  KE  is  in  the  area  to  the  west  of  Point  Sur. 
Once  (ii)  relaxation  occurs  (21  August),  energy  transfer  in  this  location  is  significantly 
reduced  and  the  location  of  energy  transfer  moves  to  the  region  north  of  Monterey  Bay, 
near  Point  Ano  Nuevo.  The  generated  mesoscale  structures  propagate  northward  along 
the  coast  with  dispersion  properties  similar  to  those  of  free  thermocline-trapped, 
coastally-trapped  waves.  These  waves  have  been  previously  observed  off  Point  Sur  by 
other  investigators  (Ramp  et  ah,  1997). 

Lagrangian  Coherent  Structures  in  the  Monterey  Bay 

Lagrangian  Coherent  Structures,  or  ECS,  is  a  computational  technique  for 
processing  velocity  field  data  in  a  flow,  including  currents  in  the  ocean  (Shadden  et  ah, 
2005).  The  basic  methodology  for  this  technique,  based  on  the  idea  of  finite  time 
Liapunov  exponents  (measuring  how  particles  are  separating  in  forward  time,  or 
converging  in  backward  time)  originated  during  the  few  years  before  AOSN-II  (Haller, 
2002;  Lekien,  2003)  and  was  subsequently  further  developed  using  the  AOSN-II  data. 
One  of  the  things  ECS  can  provide  is  to  discern  regions  of  different  particle  (or  drifter) 
behavior.  In  conjunction  with  model  predictions,  this  can  be  a  useful  predictive  tool.  The 
interesting  fact  is  that  even  in  a  complex  flow  such  as  that  in  Monterey  Bay,  ECS 
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provides  time-varying  eurves  that  separate  partieles  that  have  different  dynamieal  fates 
and  partieles  do  not  eross  these  LCS.  There  is,  for  instanee,  a  well  defined  eurve 
(teehnieally,  a  repelling  LCS)  that  separates  partieles  in  the  inner  hay  from  those  in  the 
outer  hay  (Figure  15).  Even  though  this  eurve  is  moving  in  time,  and  even  though  the 
flow  is  complex  with  multiple  space  and  time  scales,  this  curve  on  the  scale  of  the  bay 
itself  sharply  separates  particles  that  have  different  fates.  Particles  on  one  side  of  the 
curve  stay  on  that  side.  Thus,  the  fluid  to  the  right  of  the  red  curve  may  be  regarded,  for 
this  time  period,  to  be  re-circulating  in  the  bay.  That  there  even  is  such  a  well  defined 
curve  is  remarkable.  Numerical  evidence  supports  the  idea  that  this  curve  (although 
perhaps  not  LCS  for  smaller  scale  structures)  is  robust  to  data  and  model  uncertainty. 

Since  the  problem  is  time-varying,  the  recirculation  zone  does  not  persist 
indefinitely,  and  how  such  structures  appear  and  disappear  is  related  to  the  mechanisms 
of  mixing  and  transport  in  the  ocean.  LCS  can  determine  the  boundaries  of  vortex  regions 
and  elucidate  the  mechanisms  of  entrainment  and  detrainment  from  vortex  structures.  It 
does  so  with  greater  clarity  and  precision  than,  for  example,  vorticity  plots,  as  has  been 
demonstrated  in  computations  with  laboratory  vortex  rings  (Shadden  et  ah,  2006).  An 
example  for  the  ocean  has  been  computed  using  output  from  the  JPL/ROMS  model 
(Figure  16).  The  fluid  in  the  lobes,  defined  as  the  regions  initially  bounded  by  the  LCS, 
was  initially  distributed  as  in  Figure  16b.  Over  time,  the  fluid  moved  away  from  the 
repelling  LCS  (orange)  and  towards  the  attracting  LCS  (blue)  ultimately  being  entrained 
by  the  blue  vorticies  as  shown  in  Ligures  16c  and  d.  Transport  and  mixing  occurs  via  the 
motion  of  the  lobes  since  no  fluid  crosses  the  LCS.  That  drifters  in  the  green  lobes  are 
flushed  out  to  sea  while  drifters  in  the  brown  lobes  recirculate  near  the  coast  evinces  the 
property  that  LCS  separate  regions  with  different  particle  fates  and  residence  times.  This 
makes  LCS  a  very  promising  tool  for  studying  the  distribution  of  water  mass  properties, 
determining  the  fate  of  pollutants  and  oil  spills,  and  planning  how  to  seed  surface  drifters 
and  autonomous  vehicles  for  field  studies  of  ocean  dynamics. 

Ocean  Leaving  Radiance 

The  ocean-leaving  radiance  data  were  sometimes  helpful  in  delineating  the  water 
mass  boundaries  in  addition  to  showing  the  distribution  of  marine  organisms,  particularly 
during  weak  winds  when  a  thin  layer  of  warm  surface  water  masked  the  underlying 
structures.  The  hyperspectral  data  were  collected  using  a  HOBI  Labs  HydroRad 
hyperspectral  radiometer  measuring  downwelling  irradiance,  (Ed),  upwelling  radiance 
(Lu),  and  downwelling  radiance  (Ld).  Lrom  these  parameters,  the  water-leaving  radiance 
(Lw)  was  calculated  using  different  algorithms  for  a  cloud-covered  or  cloud-free  sky  as 
appropriate.  The  water-leaving  radiance  from  four  spectral  channels  was  combined 
subsequently  to  estimate  the  chlorophyll-a  (Chl-a)  content  of  the  surface  water  using  the 
SeaWiFS  OC4  version  4  algorithm.  Details  of  the  calculations  can  be  found  in  a  separate 
technical  report  (Anderson  and  Ramp,  2005). 
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Three  examples  of  the  surfaee  distributions  of  Chl-a  from  Aug.  15,  20,  and  29 
were  sampled  during  the  first  upwelling  event,  the  wind  relaxation,  and  the  seeond 
upwelling  event  (Figure  17).  Comparison  with  Figures  5g,  h,  and  1  shows  almost  no 
ehlorophyll-a  seaward  of  the  17  °C  isotherm,  whieh  may  be  taken  as  the  boundary 
between  the  California  Current  water  and  the  eoastal  water.  The  exeeption  was  an 
offshore  pateh  of  high  ehlorophyll  near  36.9°N,  123. 0°W  on  August  29  whieh  aeeording 
to  SST  (see  Figure  5k)  was  loeated  in  old  upwelled  water  whieh  had  been  previously 
adveeted  offshore  on  August  25.  All  the  Chl-a  maps  were  patchier  than  SST  with  the 
highest  values  occurring  away  from  the  coast  in  older  upwelled  water.  The  vigorous 
advection  and  mixing  during  active  upwelling  events  apparently  prohibited  the  static 
stability  required  for  a  phytoplankton  bloom  to  take  place.  Once  the  wind  died  down  and 
surface  warming  re-stratified  the  water  column,  growth  accelerated  rapidly.  The  offshore 
“squirt”  south  of  the  bay  (Figure  5h)  contained  patches  of  high  Chl-a  water  being 
adveeted  offshore  from  Point  Pinos  by  the  eddy  identified  by  the  gliders.  The  Chl-a 
distributions  on  August  15  and  August  29  show  the  greater  offshore  advection  of  material 
during  the  earlier  event  vs.  the  later  one.  The  nearly  oligotrophic  water  close  to  shore  off 
Ano  Nuevo  on  August  29  (dark  blue  in  Figure  17c)  further  supports  the  notion  that  the 
California  Current  confined  the  upwelling  center  to  the  near  shore  region  during  the 
second  event  and  encouraged  the  southward,  rather  than  offshore  advection  of  the  cold, 
nutrient-rich  water. 

Impacts  and  Conclusions 

The  AOSN-II  experiment  successfully  demonstrated  the  viability  of  several  new 
technologies  and  techniques  in  addition  to  illuminating  some  basic  scientific  issues 
surrounding  coastal  upwelling  off  central  California.  Gliders  running  closed  cells  and 
racetracks  were  able  to  determine  robust  statistics  of  the  flow  due  to  their  impressive 
persistence  and  spatial  coverage.  They  also  observed  smaller-scale  features  such  as  the 
eddy  at  the  mouth  of  the  Monterey  Bay.  Gliders  flying  formations  were  able  to  sample 
frontal  gradients  in  both  the  horizontal  and  vertical  while  moving  across  or  against  the 
flow  field  as  well  as  moving  with  it.  These  observations  were  not  previously  possible 
with  floats  or  drifters,  and  were  not  possible  without  an  impressive  level  of  autonomous 
control  over  the  small  glider  fleets.  The  low-flying  aircraft  was  able  to  map  out  the  three- 
dimensional  structure  of  the  wind  fields  with  a  level  of  detail  sufficient  to  make 
quantitative  comparisons  with  the  latest  high-resolution  atmospheric  models.  By  flying 
beneath  the  clouds,  the  aircraft  also  collected  a  time  series  of  synoptic  snapshots  of  the 
surface  temperature  and  ocean  leaving  radiance  fields  which  were  not  available  any  other 
way.  The  SST  was  a  key  parameter  for  assimilation  into  all  three  numerical  models  on  a 
daily  basis.  The  Lagrangian  Coherent  Structures  (LCS)  analysis  technique,  relatively  new 
to  the  field  of  oceanography,  showed  great  promise  for  understanding  the  flow  dynamics 
and  particle  fates,  whether  using  observational  data  (surface  current  maps  from  HF  radar) 
or  numerical  model  (JPL/ROMS)  output. 

Using  the  combined  power  of  these  new  methodologies,  the  following  picture  of 
the  circulation  in  and  around  the  Monterey  Bay  during  a  series  of  upwelling  and 
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relaxation  events  has  emerged.  Two  upwelling  events  were  observed,  one  lasting  fourteen 
days  and  the  second  nine,  with  a  three-day  relaxation  event  in  between.  The  wind  speed 
was  similar  between  the  two  upwelling  events,  blowing  steadily  from  the  northwest  at  10 
-  15  m  s'f  During  the  relaxation  winds  were  poleward  at  about  5  m  s'h  All  three  events 
were  well  sampled  by  the  aircraft,  gliders,  and  ships. 

Both  upwelling  events  began  when  a  pool  of  cold  water  less  than  13  °C  appeared 
at  the  surface  off  Cape  Ano  Nuevo.  As  the  wind  continued  to  blow  for  about  five  days, 
the  cold  pool  expanded  and  spread  offshore  and  to  the  south  across  the  mouth  of  the 
Monterey  Bay.  A  key  difference  between  the  two  events  was  the  degree  of  offshore 
spreading:  The  first  event  spread  offshore  and  to  the  south  while  the  second  event  spread 
exclusively  to  the  south.  Both  models  and  observations  indicate  that  the  key  difference 
between  these  events  was  caused  by  the  position  of  the  offshore  meanders  and  eddies  of 
the  California  Current  System.  During  the  intervening  relaxation  event,  the  boundary 
between  the  coastal  and  CCS  waters,  as  indicated  by  the  17  °C  isotherm  and  the  Chl-a 
content,  moved  40  km  onshore  to  within  5  km  of  the  coast  at  its  closed  incursion.  This 
frontal  zone  effectively  blocked  the  offshore  advection  of  cold  water  during  the  second 
upwelling  event,  and  the  equatorward  geostrophic  flow  along  the  front  aided  the 
southward  advection  of  the  cold  water  by  the  wind.  The  plume  thinned  vertically  as  it 
moved  southward  across  the  bay  to  less  than  30  m  at  the  southern  side.  The  dynamics  of 
the  eddy/frontal  movements  remains  a  mystery  at  they  were  not  always  correlated  with 
the  local  wind  stress.  While  the  front  clearly  moved  onshore  during  the  relaxation,  it  also 
moved  onshore  at  other  times  and  did  not  retreat  offshore  during  the  second  upwelling 
event.  Using  the  entire  glider  data  set,  the  temporally-  and  spatially-lagged  data  indicate 
e-folding  scales  of  2  days  and  15  km  respectively  for  these  eddies,  much  smaller  than  the 
scales  (30  days  and  100  km)  typically  observed  in  deep-water  eddies  offshore.  The 
California  Undercurrent  was  present  in  approximately  the  same  location  throughout  the 
experiment,  further  offshore  in  the  north  than  the  south.  The  changes  in  the  vertically- 
averaged  currents  were  due  primarily  to  changes  in  the  upper  water  column  rather  than 
changes  in  the  CUC. 

These  results  suggest  directions  for  future  dynamical  studies.  The  spatial 
distribution  of  the  wind  stress  curl  suggests  that  positive  curl  on  the  inshore  side  of  a 
submesocale  atmospheric  jet  may  be  important  in  driving  upwelling  in  the  Monterey  Bay, 
in  addition  to  divergence  at  the  coast.  The  HOPS  model  suggests  that  the  energy  transfer 
from  the  large-scale  to  the  mesoscale  takes  place  via  mixed  barotropic/baroclinic 
instability  off  Point  Sur  during  upwelling  events  and  in  the  northern  Monterey  Bay 
during  relaxations.  More  focused  observational  studies  are  needed  to  verify  these  ideas. 
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Figure  1.  Sampling  locations  for  autonomous  vehicles  (left)  and  other  assets  such  as 
ships,  aircraft,  and  HF  radars  (right)  during  AOSN-Il.  Green  tracks  were  obtained  from 
the  Spray  vehicles  and  the  white  tracks  from  the  Slocums.  The  blue  “V”  track  was 
occupied  at  night  by  the  CalPoly  REMUS  vehicle.  The  other  red  tracks  were  occupied  by 
the  MBARl  DORADO  vehicle  and  the  NFS  REMUS  vehicle.  The  dark  blue  track  on  the 
right  is  the  aircraft  flight  path,  and  the  thin  dotted  line  with  station  locations  is  the 
POINT  SUR  ship  survey.  The  vectors  in  the  Monterey  Bay  are  the  surface  currents  as 
sensed  by  coastal  HF  radars.  Both  sampling  schemes  are  overlaid  on  satellite  AVHRR  sea 
surface  temperature  (SST). 
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Figure  2.  Low-pass  filtered  surface  wind  vectors  from  Buoy  M2,  located  at  36°  42'  N, 
122°  23'  W,  maintained  and  operated  by  the  Monterey  Bay  Aquarium  Research  Institute 
(MBARI).  The  gray  vertical  bars  indicate  the  time  of  the  NFS  TWIN  OTTER  aircraft 
over-flights. 
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Figure  3.  A  comparison  of  the  wind,  wind  stress,  and  the  wind  stress  curl  between  the 
observed  values  sampled  by  the  aircraft  and  the  output  from  the  3-km  COAMPS^’^  model 
on  August  15,  2003  during  an  upwelling  event.  Common  scales  have  been  adopted  for 
both  pictures.  The  wind  speed  displayed  as  vectors  are  overlaid  on  the  color-coded  wind 
stress  as  indicated  by  the  left  side  of  the  color  bar.  The  wind  stress  curl  translated  into 
upwelling  velocities  in  m  d'^  are  indicated  by  the  lower  panel  dots  and  the  right-hand 
scale  of  the  color  bar.  The  black  circle  in  the  lower  panels  indicates  the  size  of  the  domain 
used  to  compute  the  curl  of  the  wind  stress. 
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Figure  4.  As  in  Figure  3  exeept  on  August  20,  2003  during  a  relaxation  event. 
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Figure  5a-f.  Time  series  of  sea  surface  temperature  (SST,  °C)  and  surface  winds  (m  s'^)  as 
sensed  from  the  aircraft  during  August  4-  13,  2003.  The  flight  path  is  indicated  by  the 
heavy  black  lines  in  panels  A  and  B  and  the  arrows  indicate  the  magnitude  and  direction 
of  the  wind  speed. 
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Figure  5g-l.  As  in  Figure  5a-f  exeept  for  August  15  to  September  4,  2003. 
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Figure  6.  Surface  temperature,  salinity,  and  current  vectors  from  the  NCOM/ICON  model 
for  August  15,  2006  (top  panel),  and  August  29,  2006  (lower  panel). 
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Figure  7.  Three  gliders  moving  northwest  in  triangular  formation  (gray  lines)  from 
1 8 ;00:00Z  August  6  into  late  morning  August  7,  2003.  Colored  dots  indieate  the 
temperature  in  °C  at  the  5  m  (top)  and  30  m  (bottom)  levels;  blaek  eireles  and  lines  show 
initial  position  and  formation;  and  open  eireles  with  a  eross  inside  triangles  illustrate  the 
path  of  the  eenter  of  the  triangle.  Red  arrows  eorrespond  to  estimates  of  down-gradient 
temperature  in  the  horizontal  plane  along  the  path  of  the  center  of  the  triangular 
formation.  The  color  contours  indicate  SST  as  sampled  by  the  TWIN  OTTER  aircraft. 
The  blue  arrows  correspond  to  currents  estimated  with  objective  analysis  using  the 
gliders'  depth  averaged  current  measurements. 
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Figure  8.  Three  gliders  moving  southeast  in  triangular  formation  (gray  lines)  during 
August  16-17,  2003.  Colored  dots  indicate  the  temperature  in  °C  at  the  5  m  (top)  and  30 
m  (bottom)  levels;  black  circles  and  lines  show  initial  position  and  formation;  and  open 
circles  with  a  cross  inside  triangles  illustrate  the  path  of  the  center  of  the  triangle.  Red 
arrows  correspond  to  estimates  of  down-gradient  temperature  in  the  horizontal  plane 
along  the  path  of  the  center  of  the  triangular  formation.  The  color  contours  indicate  SST 
as  sampled  by  the  TWIN  OTTER  aircraft.  The  blue  arrows  correspond  to  currents 
estimated  with  objective  analysis  using  the  gliders’  depth  averaged  current  measurements. 
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Figure  9.  Raw  observations  of  5-m  temperature  (eolored  squares,  seale  to  the  right)  from 
AOSN-II  ship  surveys,  moorings  and  gliders  and  depth-averaged  veloeity  from  gliders 
(black  vectors,  scale  at  upper  right)  from  a)  July  3 1  through  August  6,  2006  and  b) 

August  8-  15,  2006.  Spray  gliders  on  lines  extending  well  offshore  sample  to  the  bottom 
or  400  m  while  Slocum  gliders  on  nearshore  generally  closed  tracks  sample  to  the  bottom 
or  200  m. 
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Figure  10.  Alongshore  and  time-averaged  temperature  measurements  (°C)  averaged  in 
10-m  depth  bins  and  4  km  bins  of  offshore  distance  (shortest  distance  to  the  coast)  for  a) 
July  31  though  August  6,  2006  and  b)  August  8-15,  2006.  The  view  is  southward  with  the 
coast  on  the  left. 
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Figure  11.  Absolute  alongshore  veloeity  (cm/s)  averaged  over  a)  July  31  through  August 
6,  2006  and  b)  August  8-15,  2006.  Geostrophic  shear  (from  T  and  S  sections  like  Figure 
10)  was  referenced  by  the  alongshore-averaged  glider-measured  depth-averaged  velocity 
(from  Figure  9).  Since  this  constrains  the  vertical  integral  of  velocity,  400  m  is  not  a 
reference  level  and  the  velocities  there  are  not  weak  by  construction.  The  velocity  scale  is 
at  right  with  positive  numbers  poleward  and  negative  equatorward.  The  view  is 
southward  with  the  coast  on  the  left. 
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Figure  12.  Temporally  lagged  eorrelation  (left)  and  spatially  lagged  isotropie  correlation 
(right)  of  temperature  at  various  depths.  Temporally  lagged  result  includes  all  AOSN-11 
observations  that  are  within  5  km  of  each  other  while  spatially  lagged  correlation 
includes  all  data  whose  observation  time  differs  by  less  than  12  hours.  The  smallest  lag 
bins  do  not  include  data  correlated  with  themselves.  The  legend  at  upper  right  describes 
the  color  and  symbol  code  for  different  depths  in  meters. 
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Figure  13.  Energy  transfer  from  the  HOPS  large-seale  window  to  mesoseale  window  at 
10m.  Top  row-  available  potential  energy;  bottom  row  -  kinetie  energy. 
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Figure  14.  Schematic  of  the  major  dynamical  processes  in  the  upper  layer  near  Monterey 
Bay  during  (a)  upwelling  favorable  winds  and  (b)  relaxation.  Major  energy  transfer  from 
the  large-scale  window  to  the  mesoscale  window  takes  place  west  of  Point  Sur  during 
upwelling,  while,  during  relaxation,  major  energy  transfer  takes  place  north  of  and  in 
Monterey  Bay.  The  generated  mesoscale  events  propagate  northward  in  the  form  of 
coastally  trapped  waves. 
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Figure  15.  Four  snapshots  of  particle  dispersion  in  the  Monterey  Bay  in  the  vicinity  of  an 
LCS  (red  line)  during  upwelling  favorable  winds.  The  LCS  and  particle  dispersion  were 
computed  from  pseudo-drifters  inserted  into  the  CODAR  HF  radar  surface  current  maps 
obtained  during  the  AOSN-11  experiment.  For  additional  discussion  of  this  figure,  see 
(Shadden  et  ah,  this  issue). 
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Figure  16.  LCS  along  the  eentral  California  eoast  near  the  Monterey  Bay  as  eomputed 
from  the  JPL/ROMS  model  output.  The  LCS  are  depleted  for  (a  and  b)  July  25,  2006  at 
0900,  (e)  July  26  at  0100,  and  (d)  July  26  at  2300.  Interseetions  of  the  attracting  LCS 
(blue  curves)  and  the  repelling  LCS  (orange  curves)  define  the  boundaries  of  lobes 
colored  green  and  brown. 
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Figure  17.  Chlorophyll-a  distributions  as  sensed  from  the  hyperspeetral  sensor  on  board 
the  TWIN  OTTER  aireraft  for  a)  August  15,  b)  August  20,  and  c)  August  29,  2003. 
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